Package
library(readxl)
library(writexl)
library(dplyr)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
library(Gifi)
## Warning: package 'Gifi' was built under R version 4.2.3
library(fclust)
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ forcats 1.0.0 ✔ readr 2.1.4
## ✔ ggplot2 3.4.1 ✔ stringr 1.5.0
## ✔ lubridate 1.9.2 ✔ tibble 3.2.0
## ✔ purrr 1.0.1 ✔ tidyr 1.3.0
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ℹ Use the ]8;;http://conflicted.r-lib.org/conflicted package]8;; to force all conflicts to become errors
library(fmsb)
## Warning: package 'fmsb' was built under R version 4.2.3
library(plotly)
##
## Attaching package: 'plotly'
##
## The following object is masked from 'package:ggplot2':
##
## last_plot
##
## The following object is masked from 'package:stats':
##
## filter
##
## The following object is masked from 'package:graphics':
##
## layout
library(RColorBrewer)
library(viridis)
## Loading required package: viridisLite
Data loading
pre_df = read_excel('Datos FAMD.xlsx')
head(pre_df)
df = pre_df
colnames(df) = paste0('var', seq(ncol(df)))
df = df %>%
mutate_all(as.factor) %>%
as.data.frame()
head(df)
Categorical principal component analysis
fitord <- princals(df, ndim = 6) ## ordinal PCA
summary(fitord)
##
## Loadings (cutoff = 0.1):
## Comp1 Comp2 Comp3 Comp4 Comp5 Comp6
## var8 0.566 -0.186 -0.386 0.188 0.362
## var9 -0.841
## var10 -0.870
## var15 0.683 0.243 0.217 -0.288
## var3 0.889 -0.319 -0.248 0.131
## var5 0.894 -0.335 -0.215 0.126
## var12 -0.291 -0.328 -0.506 -0.104 0.503
## var13 -0.407 -0.375 -0.635
## var14 0.516 -0.207 -0.577 0.121
## var4 -0.178 0.294 -0.355 0.622 -0.175 -0.135
## var6 -0.251 0.128 -0.174 0.628 -0.119
## var7 -0.110 0.509 0.460
## var1 -0.210 -0.192 0.770 -0.148
## var2 0.124 -0.162 -0.108 -0.311 -0.597 -0.106
## var11 -0.275 -0.224 -0.494 -0.149 -0.662
##
## Importance (Variance Accounted For):
## Comp1 Comp2 Comp3 Comp4 Comp5 Comp6
## Eigenvalues 2.9680 2.2113 1.7880 1.4077 1.2680 0.9836
## VAF 19.7869 14.7420 11.9198 9.3848 8.4533 6.5574
## Cumulative VAF 19.7900 34.5300 46.4500 55.8300 64.2900 70.8400
# Biplot of 1 and 2 components
plot(fitord, "loadplot", plot.dim = c(1,2))
# Percentage explained variance by component
vaf = 100 * fitord$evals/sum(fitord$evals)
graf_scree = ggplot()+
aes(x=as.factor(seq(15)), y=vaf)+
geom_col(fill = '#4682b4')+
geom_text(aes(label = paste0(round(vaf,1), '%')), nudge_y = 1)+
labs(x = 'Component',
y = 'Percentage of the variance explained',
title = 'Scree Plot PRINCAL')+
theme_bw()
graf_scree
Fuzzy k-means clustering
set.seed(2022)
ran_c = seq(2, 7, 1)
ran_m = seq(1.5, 1.9, .1)
xb_ = matrix(NA, nrow = length(ran_c), ncol = length(ran_m))
for(i in ran_c){
ms = NULL
for(j in ran_m){
clust = FKM(fitord$objectscores, k = i, m = j)
# Xie and Beni index
xb = XB(clust$Xca,clust$U,clust$H,clust$m)
ms = c(ms, xb)
}
xb_[i-1,]=ms
}
dimnames(xb_) = list('cluster'=ran_c, 'fuzziness'=ran_m)
round(xb_, 3)
## fuzziness
## cluster 1.5 1.6 1.7 1.8 1.9
## 2 2.071 2.128 2.204 2.309 2.456
## 3 1.231 1.374 1.631 2.071 2.888
## 4 1.287 0.793 1.242 1.177 1.163
## 5 0.932 0.929 0.911 0.872 0.821
## 6 1.147 1.517 1.377 0.968 0.846
## 7 1.302 1.232 1.082 0.947 0.839
mej_xb = which(xb_ == min(xb_), arr.ind = TRUE)
ran_c[mej_xb[1]]
## [1] 4
ran_m[mej_xb[2]]
## [1] 1.6
plot_ly(type = 'mesh3d',
x = rep(ran_c, length(ran_m)),
y = rep(ran_m, each = length(ran_c)),
z = as.vector(xb_),
intensity = as.vector(xb_),
color = as.vector(xb_),
colors = viridis(30)) %>%
layout(scene = list(xaxis=list(title='# cluster'),
yaxis=list(title='fuzziness'),
zaxis=list(title='Xie-Beni')))
## Warning: 'mesh3d' objects don't have these attributes: 'marker'
## Valid attributes include:
## 'alphahull', 'autocolorscale', 'cauto', 'cmax', 'cmid', 'cmin', 'color', 'coloraxis', 'colorbar', 'colorscale', 'contour', 'customdata', 'customdatasrc', 'delaunayaxis', 'facecolor', 'facecolorsrc', 'flatshading', 'hoverinfo', 'hoverinfosrc', 'hoverlabel', 'hovertemplate', 'hovertemplatesrc', 'hovertext', 'hovertextsrc', 'i', 'ids', 'idssrc', 'intensity', 'intensitymode', 'intensitysrc', 'isrc', 'j', 'jsrc', 'k', 'ksrc', 'legendgroup', 'legendgrouptitle', 'legendrank', 'lighting', 'lightposition', 'meta', 'metasrc', 'name', 'opacity', 'reversescale', 'scene', 'showlegend', 'showscale', 'stream', 'text', 'textsrc', 'type', 'uid', 'uirevision', 'vertexcolor', 'vertexcolorsrc', 'visible', 'x', 'xcalendar', 'xhoverformat', 'xsrc', 'y', 'ycalendar', 'yhoverformat', 'ysrc', 'z', 'zcalendar', 'zhoverformat', 'zsrc', 'key', 'set', 'frame', 'transforms', '_isNestedKey', '_isSimpleKey', '_isGraticule', '_bbox'
# Fuzzy with best k
set.seed(123)
clust = FKM(fitord$objectscores, k = 4, m = 1.6, index = 'XB')
grupo = clust$clus[,1]
table(grupo)
## grupo
## 1 2 3 4
## 43 29 23 36
df_res = df
colnames(df_res) = colnames(pre_df)
dfsl = list(
total = cbind(df_res, grupo),
grupo1 = df_res[grupo==1,],
grupo2 = df_res[grupo==2,],
grupo3 = df_res[grupo==3,],
grupo4 = df_res[grupo==4,]
)
dfsl
## $total
## Sexo Rango/edad Sabe leer y escribir TÃtulo de Educación Superior
## 1 1 4 1 0
## 2 1 2 1 2
## 3 1 4 1 0
## 4 1 4 1 0
## 5 1 3 1 0
## 6 1 3 1 2
## 7 1 3 1 2
## 8 1 2 1 0
## 9 1 3 1 2
## 10 1 2 1 0
## 11 1 4 0 0
## 12 1 3 1 0
## 13 1 4 1 2
## 14 1 5 1 0
## 15 1 4 1 0
## 16 1 4 1 0
## 17 1 3 1 2
## 18 1 3 1 1
## 19 1 2 1 2
## 20 1 2 1 2
## 21 1 2 1 0
## 22 2 3 1 2
## 23 2 3 1 0
## 24 1 5 1 0
## 25 1 5 1 0
## 26 1 5 1 0
## 27 1 5 1 0
## 28 2 4 1 0
## 29 2 4 1 0
## 30 2 3 1 0
## 31 1 4 1 1
## 32 1 4 1 2
## 33 2 4 1 0
## 34 1 4 1 2
## 35 1 4 1 2
## 36 2 3 1 1
## 37 1 4 1 0
## 38 2 2 0 0
## 39 1 4 1 0
## 40 1 3 1 1
## 41 1 2 1 0
## 42 1 2 1 0
## 43 1 4 1 2
## 44 2 4 1 0
## 45 1 5 1 0
## 46 1 5 0 0
## 47 1 3 1 0
## 48 1 4 1 0
## 49 2 2 0 0
## 50 1 3 1 2
## 51 1 3 1 2
## 52 1 3 1 0
## 53 1 2 1 0
## 54 1 2 1 0
## 55 2 2 1 0
## 56 1 2 1 0
## 57 1 4 1 2
## 58 1 2 1 2
## 59 1 4 1 0
## 60 1 4 1 1
## 61 2 2 1 0
## 62 1 2 1 0
## 63 1 2 1 2
## 64 2 2 1 0
## 65 2 2 1 1
## 66 2 3 1 1
## 67 1 3 1 2
## 68 1 5 1 1
## 69 1 3 0 0
## 70 1 2 1 0
## 71 1 3 1 0
## 72 1 3 1 1
## 73 2 3 1 0
## 74 2 2 1 1
## 75 2 2 1 2
## 76 1 2 1 0
## 77 1 3 1 0
## 78 1 2 1 2
## 79 1 4 1 2
## 80 1 2 1 2
## 81 2 4 1 0
## 82 1 3 1 0
## 83 2 4 1 1
## 84 1 2 1 0
## 85 1 3 1 2
## 86 1 3 1 1
## 87 2 2 1 0
## 88 1 3 1 0
## 89 1 3 1 0
## 90 1 4 1 0
## 91 2 2 1 0
## 92 1 5 1 0
## 93 1 5 1 2
## 94 1 4 1 2
## 95 1 2 1 2
## 96 2 3 1 0
## 97 2 4 1 0
## 98 1 3 1 0
## 99 2 2 1 1
## 100 1 3 1 1
## 101 1 3 1 1
## 102 2 3 1 2
## 103 1 3 1 1
## 104 1 3 1 0
## 105 1 5 1 1
## 106 1 2 0 0
## 107 1 2 1 1
## 108 1 4 1 0
## 109 1 1 1 0
## 110 1 4 1 1
## 111 1 4 1 2
## 112 2 4 1 1
## 113 1 3 1 2
## 114 1 3 1 0
## 115 1 3 1 0
## 116 1 2 1 0
## 117 1 3 1 0
## 118 1 3 1 2
## 119 1 4 1 0
## 120 1 3 1 2
## 121 1 3 1 2
## 122 1 3 1 0
## 123 1 2 1 0
## 124 1 4 1 2
## 125 1 3 1 2
## 126 2 2 1 1
## 127 1 4 1 2
## 128 1 3 1 0
## 129 1 2 1 0
## 130 1 2 1 2
## 131 2 4 1 0
## Nivel de Educación Categoria de productores - Fedepalma
## 1 1 1
## 2 2 2
## 3 1 1
## 4 1 1
## 5 1 1
## 6 2 1
## 7 2 1
## 8 1 1
## 9 2 1
## 10 1 1
## 11 0 1
## 12 1 1
## 13 2 1
## 14 1 1
## 15 1 1
## 16 1 1
## 17 2 1
## 18 2 1
## 19 2 1
## 20 2 2
## 21 1 1
## 22 2 1
## 23 1 1
## 24 1 1
## 25 1 1
## 26 1 1
## 27 1 1
## 28 1 1
## 29 1 1
## 30 1 1
## 31 2 1
## 32 2 1
## 33 1 1
## 34 2 1
## 35 2 1
## 36 2 1
## 37 1 1
## 38 0 1
## 39 1 1
## 40 2 1
## 41 1 1
## 42 1 1
## 43 2 1
## 44 1 1
## 45 1 1
## 46 0 1
## 47 1 1
## 48 1 1
## 49 0 1
## 50 2 1
## 51 2 1
## 52 1 1
## 53 1 1
## 54 1 1
## 55 1 1
## 56 1 1
## 57 2 1
## 58 2 1
## 59 1 1
## 60 2 1
## 61 1 1
## 62 1 1
## 63 2 1
## 64 1 1
## 65 2 1
## 66 2 1
## 67 2 2
## 68 2 1
## 69 0 1
## 70 1 1
## 71 1 1
## 72 2 1
## 73 1 1
## 74 2 1
## 75 2 1
## 76 1 1
## 77 1 1
## 78 2 1
## 79 2 2
## 80 2 2
## 81 1 1
## 82 1 1
## 83 2 1
## 84 1 1
## 85 2 1
## 86 2 1
## 87 1 1
## 88 1 1
## 89 1 1
## 90 1 1
## 91 1 1
## 92 1 1
## 93 2 1
## 94 2 1
## 95 2 1
## 96 1 1
## 97 1 1
## 98 1 1
## 99 2 1
## 100 2 1
## 101 2 1
## 102 2 1
## 103 2 1
## 104 1 1
## 105 2 1
## 106 0 1
## 107 2 1
## 108 1 1
## 109 1 1
## 110 2 1
## 111 2 1
## 112 2 1
## 113 2 1
## 114 1 1
## 115 1 1
## 116 1 1
## 117 1 1
## 118 2 1
## 119 1 1
## 120 2 1
## 121 2 1
## 122 1 1
## 123 1 1
## 124 2 1
## 125 2 1
## 126 2 1
## 127 2 1
## 128 1 1
## 129 1 1
## 130 2 1
## 131 1 1
## Adopción de tecnologÃa dicotómica
## 1 1
## 2 1
## 3 2
## 4 1
## 5 1
## 6 2
## 7 1
## 8 1
## 9 1
## 10 1
## 11 1
## 12 1
## 13 1
## 14 1
## 15 1
## 16 2
## 17 1
## 18 1
## 19 1
## 20 1
## 21 1
## 22 1
## 23 2
## 24 1
## 25 1
## 26 1
## 27 1
## 28 2
## 29 2
## 30 1
## 31 1
## 32 1
## 33 1
## 34 1
## 35 2
## 36 1
## 37 1
## 38 2
## 39 1
## 40 2
## 41 1
## 42 1
## 43 1
## 44 1
## 45 1
## 46 1
## 47 1
## 48 1
## 49 1
## 50 1
## 51 1
## 52 1
## 53 1
## 54 1
## 55 1
## 56 1
## 57 2
## 58 1
## 59 1
## 60 1
## 61 1
## 62 1
## 63 1
## 64 1
## 65 1
## 66 2
## 67 2
## 68 1
## 69 1
## 70 1
## 71 1
## 72 1
## 73 1
## 74 2
## 75 2
## 76 1
## 77 1
## 78 2
## 79 2
## 80 1
## 81 1
## 82 1
## 83 1
## 84 1
## 85 1
## 86 1
## 87 1
## 88 1
## 89 1
## 90 1
## 91 1
## 92 1
## 93 1
## 94 1
## 95 1
## 96 1
## 97 1
## 98 1
## 99 1
## 100 2
## 101 1
## 102 1
## 103 1
## 104 1
## 105 1
## 106 2
## 107 2
## 108 1
## 109 1
## 110 1
## 111 1
## 112 1
## 113 1
## 114 2
## 115 1
## 116 1
## 117 1
## 118 1
## 119 1
## 120 1
## 121 1
## 122 1
## 123 1
## 124 1
## 125 1
## 126 1
## 127 1
## 128 1
## 129 1
## 130 1
## 131 1
## Trabaja en actividades diferentes al cultivo de palma Tipo de tenencia
## 1 NO 1
## 2 SI 2
## 3 NO 3
## 4 SI 2
## 5 NO 3
## 6 NO 2
## 7 SI 2
## 8 SI 3
## 9 NO 3
## 10 SI 2
## 11 NO 3
## 12 SI 2
## 13 NO 2
## 14 NO 2
## 15 SI 3
## 16 NO 1
## 17 NO 3
## 18 NO 2
## 19 NO 3
## 20 NO 3
## 21 SI 2
## 22 NO 3
## 23 NO 3
## 24 SI 2
## 25 NO 2
## 26 NO 2
## 27 NO 2
## 28 NO 2
## 29 NO 2
## 30 NO 1
## 31 SI 1
## 32 SI 2
## 33 NO 2
## 34 NO 3
## 35 NO 3
## 36 SI 1
## 37 NO 3
## 38 NO 3
## 39 NO 1
## 40 SI 2
## 41 NO 2
## 42 NO 3
## 43 SI 2
## 44 NO 2
## 45 SI 2
## 46 SI 2
## 47 NO 2
## 48 NO 2
## 49 SI 2
## 50 NO 2
## 51 SI 2
## 52 NO 2
## 53 NO 4
## 54 NO 2
## 55 NO 2
## 56 NO 2
## 57 SI 2
## 58 SI 2
## 59 SI 2
## 60 SI 2
## 61 NO 2
## 62 NO 2
## 63 NO 3
## 64 NO 3
## 65 NO 3
## 66 NO 3
## 67 NO 3
## 68 NO 3
## 69 SI 2
## 70 SI 2
## 71 NO 2
## 72 NO 3
## 73 NO 3
## 74 SI 3
## 75 SI 2
## 76 NO 2
## 77 SI 2
## 78 SI 2
## 79 NO 3
## 80 NO 3
## 81 SI 2
## 82 NO 3
## 83 NO 2
## 84 SI 2
## 85 NO 2
## 86 NO 2
## 87 NO 1
## 88 NO 2
## 89 NO 2
## 90 NO 1
## 91 NO 2
## 92 NO 3
## 93 NO 3
## 94 NO 2
## 95 SI 2
## 96 SI 2
## 97 SI 2
## 98 NO 2
## 99 NO 2
## 100 NO 2
## 101 NO 3
## 102 NO 3
## 103 NO 3
## 104 NO 2
## 105 SI 2
## 106 NO 2
## 107 NO 3
## 108 NO 2
## 109 NO 3
## 110 NO 3
## 111 NO 3
## 112 NO 2
## 113 NO 3
## 114 NO 2
## 115 NO 3
## 116 NO 3
## 117 NO 3
## 118 SI 2
## 119 NO 3
## 120 NO 3
## 121 NO 3
## 122 NO 3
## 123 NO 1
## 124 NO 2
## 125 NO 3
## 126 NO 3
## 127 NO 3
## 128 NO 3
## 129 NO 3
## 130 SI 2
## 131 SI 5
## Uso del predio Utiliza sistema de riego Utiliza fertilizantes
## 1 1 SI SI
## 2 1 SI SI
## 3 2 SI SI
## 4 2 SI SI
## 5 2 SI SI
## 6 2 SI SI
## 7 1 SI SI
## 8 1 SI SI
## 9 2 SI SI
## 10 1 SI SI
## 11 2 SI SI
## 12 2 SI SI
## 13 1 SI SI
## 14 2 SI SI
## 15 2 SI SI
## 16 1 NO SI
## 17 2 SI SI
## 18 1 SI NO
## 19 2 SI SI
## 20 2 SI SI
## 21 2 SI SI
## 22 2 SI SI
## 23 2 SI SI
## 24 1 SI SI
## 25 1 SI SI
## 26 1 SI SI
## 27 1 SI SI
## 28 1 SI SI
## 29 1 SI SI
## 30 1 SI SI
## 31 1 SI SI
## 32 1 SI SI
## 33 2 SI SI
## 34 2 SI SI
## 35 2 SI SI
## 36 1 SI SI
## 37 2 SI SI
## 38 2 SI SI
## 39 2 SI SI
## 40 1 SI SI
## 41 1 SI SI
## 42 2 SI SI
## 43 1 SI SI
## 44 1 SI NO
## 45 1 SI SI
## 46 1 SI SI
## 47 1 SI SI
## 48 1 SI SI
## 49 1 SI SI
## 50 1 SI NO
## 51 1 SI SI
## 52 1 SI NO
## 53 1 SI SI
## 54 1 SI SI
## 55 1 SI SI
## 56 1 SI SI
## 57 1 SI SI
## 58 1 SI SI
## 59 1 SI SI
## 60 1 SI SI
## 61 1 SI SI
## 62 1 SI SI
## 63 2 SI SI
## 64 2 SI SI
## 65 2 SI SI
## 66 2 SI SI
## 67 2 SI SI
## 68 2 SI SI
## 69 1 SI SI
## 70 1 NO NO
## 71 1 SI SI
## 72 2 SI SI
## 73 2 SI SI
## 74 1 SI SI
## 75 1 SI SI
## 76 1 SI SI
## 77 1 NO SI
## 78 1 SI SI
## 79 2 SI SI
## 80 2 SI SI
## 81 1 SI SI
## 82 2 SI SI
## 83 1 SI SI
## 84 1 SI SI
## 85 1 SI SI
## 86 1 SI SI
## 87 1 SI SI
## 88 1 SI SI
## 89 1 SI SI
## 90 1 SI SI
## 91 1 SI SI
## 92 2 SI SI
## 93 1 SI NO
## 94 1 SI SI
## 95 1 SI SI
## 96 1 SI SI
## 97 1 SI SI
## 98 1 SI SI
## 99 2 SI SI
## 100 1 SI NO
## 101 2 SI SI
## 102 2 SI SI
## 103 2 SI SI
## 104 2 SI SI
## 105 1 SI SI
## 106 1 SI SI
## 107 2 SI SI
## 108 1 SI SI
## 109 2 SI SI
## 110 2 SI SI
## 111 2 SI SI
## 112 1 SI SI
## 113 2 SI SI
## 114 2 SI NO
## 115 2 SI SI
## 116 2 SI SI
## 117 2 SI SI
## 118 1 SI SI
## 119 2 SI SI
## 120 2 SI SI
## 121 2 SI SI
## 122 2 SI SI
## 123 1 SI NO
## 124 1 SI SI
## 125 2 SI SI
## 126 2 SI SI
## 127 2 SI SI
## 128 2 SI SI
## 129 2 SI SI
## 130 1 SI SI
## 131 1 SI SI
## Posee esquema de manejo de plagas y enfermedades
## 1 NO
## 2 SI
## 3 SI
## 4 SI
## 5 SI
## 6 SI
## 7 SI
## 8 SI
## 9 SI
## 10 SI
## 11 SI
## 12 SI
## 13 SI
## 14 SI
## 15 SI
## 16 NO
## 17 SI
## 18 NO
## 19 SI
## 20 SI
## 21 SI
## 22 SI
## 23 SI
## 24 SI
## 25 SI
## 26 SI
## 27 SI
## 28 SI
## 29 SI
## 30 NO
## 31 SI
## 32 NO
## 33 SI
## 34 SI
## 35 SI
## 36 SI
## 37 SI
## 38 SI
## 39 SI
## 40 SI
## 41 SI
## 42 SI
## 43 SI
## 44 NO
## 45 SI
## 46 SI
## 47 SI
## 48 SI
## 49 SI
## 50 SI
## 51 SI
## 52 NO
## 53 SI
## 54 SI
## 55 SI
## 56 SI
## 57 SI
## 58 SI
## 59 SI
## 60 SI
## 61 SI
## 62 SI
## 63 SI
## 64 SI
## 65 SI
## 66 SI
## 67 SI
## 68 SI
## 69 SI
## 70 NO
## 71 SI
## 72 SI
## 73 SI
## 74 SI
## 75 SI
## 76 NO
## 77 NO
## 78 SI
## 79 SI
## 80 SI
## 81 SI
## 82 SI
## 83 SI
## 84 SI
## 85 NO
## 86 SI
## 87 SI
## 88 NO
## 89 SI
## 90 NO
## 91 SI
## 92 SI
## 93 SI
## 94 SI
## 95 SI
## 96 SI
## 97 SI
## 98 SI
## 99 SI
## 100 SI
## 101 SI
## 102 SI
## 103 SI
## 104 SI
## 105 SI
## 106 SI
## 107 SI
## 108 SI
## 109 SI
## 110 SI
## 111 SI
## 112 SI
## 113 SI
## 114 NO
## 115 SI
## 116 SI
## 117 SI
## 118 SI
## 119 SI
## 120 SI
## 121 SI
## 122 SI
## 123 SI
## 124 SI
## 125 SI
## 126 SI
## 127 SI
## 128 SI
## 129 SI
## 130 SI
## 131 SI
## Ganancias del cultivo de palma Tipo de trabajadores grupo
## 1 2 1 1
## 2 4 2 3
## 3 2 1 4
## 4 3 1 1
## 5 2 1 4
## 6 3 1 3
## 7 1 2 3
## 8 3 2 3
## 9 2 1 4
## 10 4 2 3
## 11 2 1 4
## 12 3 1 1
## 13 2 2 3
## 14 3 1 1
## 15 4 1 1
## 16 1 1 1
## 17 2 1 4
## 18 1 3 1
## 19 2 1 4
## 20 2 1 4
## 21 3 1 1
## 22 2 1 4
## 23 2 1 2
## 24 3 1 1
## 25 3 2 1
## 26 2 1 1
## 27 5 1 1
## 28 3 2 2
## 29 3 2 2
## 30 1 2 2
## 31 2 3 3
## 32 3 2 3
## 33 4 1 2
## 34 2 1 4
## 35 2 1 4
## 36 1 1 2
## 37 2 1 4
## 38 2 1 2
## 39 2 3 1
## 40 3 1 3
## 41 2 1 1
## 42 2 1 4
## 43 2 2 3
## 44 1 2 2
## 45 4 2 1
## 46 1 2 1
## 47 3 1 1
## 48 3 1 1
## 49 5 1 2
## 50 1 2 3
## 51 4 1 3
## 52 1 2 1
## 53 3 2 1
## 54 3 2 1
## 55 3 1 2
## 56 2 1 1
## 57 5 2 3
## 58 5 2 3
## 59 2 1 1
## 60 3 1 1
## 61 4 1 2
## 62 2 1 1
## 63 2 1 4
## 64 2 1 2
## 65 2 1 2
## 66 2 1 2
## 67 2 1 3
## 68 2 1 4
## 69 1 1 1
## 70 2 3 1
## 71 2 1 1
## 72 2 1 4
## 73 2 1 2
## 74 3 2 2
## 75 2 2 2
## 76 2 3 1
## 77 1 3 1
## 78 5 2 3
## 79 2 1 3
## 80 2 1 4
## 81 1 1 2
## 82 2 1 4
## 83 2 1 2
## 84 3 1 1
## 85 3 2 3
## 86 2 1 1
## 87 3 1 2
## 88 2 1 1
## 89 3 2 1
## 90 1 3 1
## 91 3 2 2
## 92 2 1 4
## 93 1 2 1
## 94 4 2 3
## 95 3 2 3
## 96 3 3 2
## 97 3 3 2
## 98 2 1 1
## 99 3 1 2
## 100 3 1 2
## 101 2 1 4
## 102 2 1 4
## 103 2 1 4
## 104 2 2 1
## 105 3 1 1
## 106 3 1 2
## 107 2 1 4
## 108 4 2 1
## 109 2 1 4
## 110 3 2 1
## 111 2 1 4
## 112 2 1 2
## 113 2 1 4
## 114 1 2 1
## 115 2 1 4
## 116 2 1 4
## 117 2 1 4
## 118 2 2 3
## 119 2 1 4
## 120 2 1 4
## 121 2 1 4
## 122 2 1 4
## 123 2 3 1
## 124 2 2 3
## 125 2 1 4
## 126 2 1 2
## 127 2 1 4
## 128 2 1 4
## 129 2 1 4
## 130 5 2 3
## 131 4 2 2
##
## $grupo1
## Sexo Rango/edad Sabe leer y escribir TÃtulo de Educación Superior
## 1 1 4 1 0
## 4 1 4 1 0
## 12 1 3 1 0
## 14 1 5 1 0
## 15 1 4 1 0
## 16 1 4 1 0
## 18 1 3 1 1
## 21 1 2 1 0
## 24 1 5 1 0
## 25 1 5 1 0
## 26 1 5 1 0
## 27 1 5 1 0
## 39 1 4 1 0
## 41 1 2 1 0
## 45 1 5 1 0
## 46 1 5 0 0
## 47 1 3 1 0
## 48 1 4 1 0
## 52 1 3 1 0
## 53 1 2 1 0
## 54 1 2 1 0
## 56 1 2 1 0
## 59 1 4 1 0
## 60 1 4 1 1
## 62 1 2 1 0
## 69 1 3 0 0
## 70 1 2 1 0
## 71 1 3 1 0
## 76 1 2 1 0
## 77 1 3 1 0
## 84 1 2 1 0
## 86 1 3 1 1
## 88 1 3 1 0
## 89 1 3 1 0
## 90 1 4 1 0
## 93 1 5 1 2
## 98 1 3 1 0
## 104 1 3 1 0
## 105 1 5 1 1
## 108 1 4 1 0
## 110 1 4 1 1
## 114 1 3 1 0
## 123 1 2 1 0
## Nivel de Educación Categoria de productores - Fedepalma
## 1 1 1
## 4 1 1
## 12 1 1
## 14 1 1
## 15 1 1
## 16 1 1
## 18 2 1
## 21 1 1
## 24 1 1
## 25 1 1
## 26 1 1
## 27 1 1
## 39 1 1
## 41 1 1
## 45 1 1
## 46 0 1
## 47 1 1
## 48 1 1
## 52 1 1
## 53 1 1
## 54 1 1
## 56 1 1
## 59 1 1
## 60 2 1
## 62 1 1
## 69 0 1
## 70 1 1
## 71 1 1
## 76 1 1
## 77 1 1
## 84 1 1
## 86 2 1
## 88 1 1
## 89 1 1
## 90 1 1
## 93 2 1
## 98 1 1
## 104 1 1
## 105 2 1
## 108 1 1
## 110 2 1
## 114 1 1
## 123 1 1
## Adopción de tecnologÃa dicotómica
## 1 1
## 4 1
## 12 1
## 14 1
## 15 1
## 16 2
## 18 1
## 21 1
## 24 1
## 25 1
## 26 1
## 27 1
## 39 1
## 41 1
## 45 1
## 46 1
## 47 1
## 48 1
## 52 1
## 53 1
## 54 1
## 56 1
## 59 1
## 60 1
## 62 1
## 69 1
## 70 1
## 71 1
## 76 1
## 77 1
## 84 1
## 86 1
## 88 1
## 89 1
## 90 1
## 93 1
## 98 1
## 104 1
## 105 1
## 108 1
## 110 1
## 114 2
## 123 1
## Trabaja en actividades diferentes al cultivo de palma Tipo de tenencia
## 1 NO 1
## 4 SI 2
## 12 SI 2
## 14 NO 2
## 15 SI 3
## 16 NO 1
## 18 NO 2
## 21 SI 2
## 24 SI 2
## 25 NO 2
## 26 NO 2
## 27 NO 2
## 39 NO 1
## 41 NO 2
## 45 SI 2
## 46 SI 2
## 47 NO 2
## 48 NO 2
## 52 NO 2
## 53 NO 4
## 54 NO 2
## 56 NO 2
## 59 SI 2
## 60 SI 2
## 62 NO 2
## 69 SI 2
## 70 SI 2
## 71 NO 2
## 76 NO 2
## 77 SI 2
## 84 SI 2
## 86 NO 2
## 88 NO 2
## 89 NO 2
## 90 NO 1
## 93 NO 3
## 98 NO 2
## 104 NO 2
## 105 SI 2
## 108 NO 2
## 110 NO 3
## 114 NO 2
## 123 NO 1
## Uso del predio Utiliza sistema de riego Utiliza fertilizantes
## 1 1 SI SI
## 4 2 SI SI
## 12 2 SI SI
## 14 2 SI SI
## 15 2 SI SI
## 16 1 NO SI
## 18 1 SI NO
## 21 2 SI SI
## 24 1 SI SI
## 25 1 SI SI
## 26 1 SI SI
## 27 1 SI SI
## 39 2 SI SI
## 41 1 SI SI
## 45 1 SI SI
## 46 1 SI SI
## 47 1 SI SI
## 48 1 SI SI
## 52 1 SI NO
## 53 1 SI SI
## 54 1 SI SI
## 56 1 SI SI
## 59 1 SI SI
## 60 1 SI SI
## 62 1 SI SI
## 69 1 SI SI
## 70 1 NO NO
## 71 1 SI SI
## 76 1 SI SI
## 77 1 NO SI
## 84 1 SI SI
## 86 1 SI SI
## 88 1 SI SI
## 89 1 SI SI
## 90 1 SI SI
## 93 1 SI NO
## 98 1 SI SI
## 104 2 SI SI
## 105 1 SI SI
## 108 1 SI SI
## 110 2 SI SI
## 114 2 SI NO
## 123 1 SI NO
## Posee esquema de manejo de plagas y enfermedades
## 1 NO
## 4 SI
## 12 SI
## 14 SI
## 15 SI
## 16 NO
## 18 NO
## 21 SI
## 24 SI
## 25 SI
## 26 SI
## 27 SI
## 39 SI
## 41 SI
## 45 SI
## 46 SI
## 47 SI
## 48 SI
## 52 NO
## 53 SI
## 54 SI
## 56 SI
## 59 SI
## 60 SI
## 62 SI
## 69 SI
## 70 NO
## 71 SI
## 76 NO
## 77 NO
## 84 SI
## 86 SI
## 88 NO
## 89 SI
## 90 NO
## 93 SI
## 98 SI
## 104 SI
## 105 SI
## 108 SI
## 110 SI
## 114 NO
## 123 SI
## Ganancias del cultivo de palma Tipo de trabajadores
## 1 2 1
## 4 3 1
## 12 3 1
## 14 3 1
## 15 4 1
## 16 1 1
## 18 1 3
## 21 3 1
## 24 3 1
## 25 3 2
## 26 2 1
## 27 5 1
## 39 2 3
## 41 2 1
## 45 4 2
## 46 1 2
## 47 3 1
## 48 3 1
## 52 1 2
## 53 3 2
## 54 3 2
## 56 2 1
## 59 2 1
## 60 3 1
## 62 2 1
## 69 1 1
## 70 2 3
## 71 2 1
## 76 2 3
## 77 1 3
## 84 3 1
## 86 2 1
## 88 2 1
## 89 3 2
## 90 1 3
## 93 1 2
## 98 2 1
## 104 2 2
## 105 3 1
## 108 4 2
## 110 3 2
## 114 1 2
## 123 2 3
##
## $grupo2
## Sexo Rango/edad Sabe leer y escribir TÃtulo de Educación Superior
## 23 2 3 1 0
## 28 2 4 1 0
## 29 2 4 1 0
## 30 2 3 1 0
## 33 2 4 1 0
## 36 2 3 1 1
## 38 2 2 0 0
## 44 2 4 1 0
## 49 2 2 0 0
## 55 2 2 1 0
## 61 2 2 1 0
## 64 2 2 1 0
## 65 2 2 1 1
## 66 2 3 1 1
## 73 2 3 1 0
## 74 2 2 1 1
## 75 2 2 1 2
## 81 2 4 1 0
## 83 2 4 1 1
## 87 2 2 1 0
## 91 2 2 1 0
## 96 2 3 1 0
## 97 2 4 1 0
## 99 2 2 1 1
## 100 1 3 1 1
## 106 1 2 0 0
## 112 2 4 1 1
## 126 2 2 1 1
## 131 2 4 1 0
## Nivel de Educación Categoria de productores - Fedepalma
## 23 1 1
## 28 1 1
## 29 1 1
## 30 1 1
## 33 1 1
## 36 2 1
## 38 0 1
## 44 1 1
## 49 0 1
## 55 1 1
## 61 1 1
## 64 1 1
## 65 2 1
## 66 2 1
## 73 1 1
## 74 2 1
## 75 2 1
## 81 1 1
## 83 2 1
## 87 1 1
## 91 1 1
## 96 1 1
## 97 1 1
## 99 2 1
## 100 2 1
## 106 0 1
## 112 2 1
## 126 2 1
## 131 1 1
## Adopción de tecnologÃa dicotómica
## 23 2
## 28 2
## 29 2
## 30 1
## 33 1
## 36 1
## 38 2
## 44 1
## 49 1
## 55 1
## 61 1
## 64 1
## 65 1
## 66 2
## 73 1
## 74 2
## 75 2
## 81 1
## 83 1
## 87 1
## 91 1
## 96 1
## 97 1
## 99 1
## 100 2
## 106 2
## 112 1
## 126 1
## 131 1
## Trabaja en actividades diferentes al cultivo de palma Tipo de tenencia
## 23 NO 3
## 28 NO 2
## 29 NO 2
## 30 NO 1
## 33 NO 2
## 36 SI 1
## 38 NO 3
## 44 NO 2
## 49 SI 2
## 55 NO 2
## 61 NO 2
## 64 NO 3
## 65 NO 3
## 66 NO 3
## 73 NO 3
## 74 SI 3
## 75 SI 2
## 81 SI 2
## 83 NO 2
## 87 NO 1
## 91 NO 2
## 96 SI 2
## 97 SI 2
## 99 NO 2
## 100 NO 2
## 106 NO 2
## 112 NO 2
## 126 NO 3
## 131 SI 5
## Uso del predio Utiliza sistema de riego Utiliza fertilizantes
## 23 2 SI SI
## 28 1 SI SI
## 29 1 SI SI
## 30 1 SI SI
## 33 2 SI SI
## 36 1 SI SI
## 38 2 SI SI
## 44 1 SI NO
## 49 1 SI SI
## 55 1 SI SI
## 61 1 SI SI
## 64 2 SI SI
## 65 2 SI SI
## 66 2 SI SI
## 73 2 SI SI
## 74 1 SI SI
## 75 1 SI SI
## 81 1 SI SI
## 83 1 SI SI
## 87 1 SI SI
## 91 1 SI SI
## 96 1 SI SI
## 97 1 SI SI
## 99 2 SI SI
## 100 1 SI NO
## 106 1 SI SI
## 112 1 SI SI
## 126 2 SI SI
## 131 1 SI SI
## Posee esquema de manejo de plagas y enfermedades
## 23 SI
## 28 SI
## 29 SI
## 30 NO
## 33 SI
## 36 SI
## 38 SI
## 44 NO
## 49 SI
## 55 SI
## 61 SI
## 64 SI
## 65 SI
## 66 SI
## 73 SI
## 74 SI
## 75 SI
## 81 SI
## 83 SI
## 87 SI
## 91 SI
## 96 SI
## 97 SI
## 99 SI
## 100 SI
## 106 SI
## 112 SI
## 126 SI
## 131 SI
## Ganancias del cultivo de palma Tipo de trabajadores
## 23 2 1
## 28 3 2
## 29 3 2
## 30 1 2
## 33 4 1
## 36 1 1
## 38 2 1
## 44 1 2
## 49 5 1
## 55 3 1
## 61 4 1
## 64 2 1
## 65 2 1
## 66 2 1
## 73 2 1
## 74 3 2
## 75 2 2
## 81 1 1
## 83 2 1
## 87 3 1
## 91 3 2
## 96 3 3
## 97 3 3
## 99 3 1
## 100 3 1
## 106 3 1
## 112 2 1
## 126 2 1
## 131 4 2
##
## $grupo3
## Sexo Rango/edad Sabe leer y escribir TÃtulo de Educación Superior
## 2 1 2 1 2
## 6 1 3 1 2
## 7 1 3 1 2
## 8 1 2 1 0
## 10 1 2 1 0
## 13 1 4 1 2
## 31 1 4 1 1
## 32 1 4 1 2
## 40 1 3 1 1
## 43 1 4 1 2
## 50 1 3 1 2
## 51 1 3 1 2
## 57 1 4 1 2
## 58 1 2 1 2
## 67 1 3 1 2
## 78 1 2 1 2
## 79 1 4 1 2
## 85 1 3 1 2
## 94 1 4 1 2
## 95 1 2 1 2
## 118 1 3 1 2
## 124 1 4 1 2
## 130 1 2 1 2
## Nivel de Educación Categoria de productores - Fedepalma
## 2 2 2
## 6 2 1
## 7 2 1
## 8 1 1
## 10 1 1
## 13 2 1
## 31 2 1
## 32 2 1
## 40 2 1
## 43 2 1
## 50 2 1
## 51 2 1
## 57 2 1
## 58 2 1
## 67 2 2
## 78 2 1
## 79 2 2
## 85 2 1
## 94 2 1
## 95 2 1
## 118 2 1
## 124 2 1
## 130 2 1
## Adopción de tecnologÃa dicotómica
## 2 1
## 6 2
## 7 1
## 8 1
## 10 1
## 13 1
## 31 1
## 32 1
## 40 2
## 43 1
## 50 1
## 51 1
## 57 2
## 58 1
## 67 2
## 78 2
## 79 2
## 85 1
## 94 1
## 95 1
## 118 1
## 124 1
## 130 1
## Trabaja en actividades diferentes al cultivo de palma Tipo de tenencia
## 2 SI 2
## 6 NO 2
## 7 SI 2
## 8 SI 3
## 10 SI 2
## 13 NO 2
## 31 SI 1
## 32 SI 2
## 40 SI 2
## 43 SI 2
## 50 NO 2
## 51 SI 2
## 57 SI 2
## 58 SI 2
## 67 NO 3
## 78 SI 2
## 79 NO 3
## 85 NO 2
## 94 NO 2
## 95 SI 2
## 118 SI 2
## 124 NO 2
## 130 SI 2
## Uso del predio Utiliza sistema de riego Utiliza fertilizantes
## 2 1 SI SI
## 6 2 SI SI
## 7 1 SI SI
## 8 1 SI SI
## 10 1 SI SI
## 13 1 SI SI
## 31 1 SI SI
## 32 1 SI SI
## 40 1 SI SI
## 43 1 SI SI
## 50 1 SI NO
## 51 1 SI SI
## 57 1 SI SI
## 58 1 SI SI
## 67 2 SI SI
## 78 1 SI SI
## 79 2 SI SI
## 85 1 SI SI
## 94 1 SI SI
## 95 1 SI SI
## 118 1 SI SI
## 124 1 SI SI
## 130 1 SI SI
## Posee esquema de manejo de plagas y enfermedades
## 2 SI
## 6 SI
## 7 SI
## 8 SI
## 10 SI
## 13 SI
## 31 SI
## 32 NO
## 40 SI
## 43 SI
## 50 SI
## 51 SI
## 57 SI
## 58 SI
## 67 SI
## 78 SI
## 79 SI
## 85 NO
## 94 SI
## 95 SI
## 118 SI
## 124 SI
## 130 SI
## Ganancias del cultivo de palma Tipo de trabajadores
## 2 4 2
## 6 3 1
## 7 1 2
## 8 3 2
## 10 4 2
## 13 2 2
## 31 2 3
## 32 3 2
## 40 3 1
## 43 2 2
## 50 1 2
## 51 4 1
## 57 5 2
## 58 5 2
## 67 2 1
## 78 5 2
## 79 2 1
## 85 3 2
## 94 4 2
## 95 3 2
## 118 2 2
## 124 2 2
## 130 5 2
##
## $grupo4
## Sexo Rango/edad Sabe leer y escribir TÃtulo de Educación Superior
## 3 1 4 1 0
## 5 1 3 1 0
## 9 1 3 1 2
## 11 1 4 0 0
## 17 1 3 1 2
## 19 1 2 1 2
## 20 1 2 1 2
## 22 2 3 1 2
## 34 1 4 1 2
## 35 1 4 1 2
## 37 1 4 1 0
## 42 1 2 1 0
## 63 1 2 1 2
## 68 1 5 1 1
## 72 1 3 1 1
## 80 1 2 1 2
## 82 1 3 1 0
## 92 1 5 1 0
## 101 1 3 1 1
## 102 2 3 1 2
## 103 1 3 1 1
## 107 1 2 1 1
## 109 1 1 1 0
## 111 1 4 1 2
## 113 1 3 1 2
## 115 1 3 1 0
## 116 1 2 1 0
## 117 1 3 1 0
## 119 1 4 1 0
## 120 1 3 1 2
## 121 1 3 1 2
## 122 1 3 1 0
## 125 1 3 1 2
## 127 1 4 1 2
## 128 1 3 1 0
## 129 1 2 1 0
## Nivel de Educación Categoria de productores - Fedepalma
## 3 1 1
## 5 1 1
## 9 2 1
## 11 0 1
## 17 2 1
## 19 2 1
## 20 2 2
## 22 2 1
## 34 2 1
## 35 2 1
## 37 1 1
## 42 1 1
## 63 2 1
## 68 2 1
## 72 2 1
## 80 2 2
## 82 1 1
## 92 1 1
## 101 2 1
## 102 2 1
## 103 2 1
## 107 2 1
## 109 1 1
## 111 2 1
## 113 2 1
## 115 1 1
## 116 1 1
## 117 1 1
## 119 1 1
## 120 2 1
## 121 2 1
## 122 1 1
## 125 2 1
## 127 2 1
## 128 1 1
## 129 1 1
## Adopción de tecnologÃa dicotómica
## 3 2
## 5 1
## 9 1
## 11 1
## 17 1
## 19 1
## 20 1
## 22 1
## 34 1
## 35 2
## 37 1
## 42 1
## 63 1
## 68 1
## 72 1
## 80 1
## 82 1
## 92 1
## 101 1
## 102 1
## 103 1
## 107 2
## 109 1
## 111 1
## 113 1
## 115 1
## 116 1
## 117 1
## 119 1
## 120 1
## 121 1
## 122 1
## 125 1
## 127 1
## 128 1
## 129 1
## Trabaja en actividades diferentes al cultivo de palma Tipo de tenencia
## 3 NO 3
## 5 NO 3
## 9 NO 3
## 11 NO 3
## 17 NO 3
## 19 NO 3
## 20 NO 3
## 22 NO 3
## 34 NO 3
## 35 NO 3
## 37 NO 3
## 42 NO 3
## 63 NO 3
## 68 NO 3
## 72 NO 3
## 80 NO 3
## 82 NO 3
## 92 NO 3
## 101 NO 3
## 102 NO 3
## 103 NO 3
## 107 NO 3
## 109 NO 3
## 111 NO 3
## 113 NO 3
## 115 NO 3
## 116 NO 3
## 117 NO 3
## 119 NO 3
## 120 NO 3
## 121 NO 3
## 122 NO 3
## 125 NO 3
## 127 NO 3
## 128 NO 3
## 129 NO 3
## Uso del predio Utiliza sistema de riego Utiliza fertilizantes
## 3 2 SI SI
## 5 2 SI SI
## 9 2 SI SI
## 11 2 SI SI
## 17 2 SI SI
## 19 2 SI SI
## 20 2 SI SI
## 22 2 SI SI
## 34 2 SI SI
## 35 2 SI SI
## 37 2 SI SI
## 42 2 SI SI
## 63 2 SI SI
## 68 2 SI SI
## 72 2 SI SI
## 80 2 SI SI
## 82 2 SI SI
## 92 2 SI SI
## 101 2 SI SI
## 102 2 SI SI
## 103 2 SI SI
## 107 2 SI SI
## 109 2 SI SI
## 111 2 SI SI
## 113 2 SI SI
## 115 2 SI SI
## 116 2 SI SI
## 117 2 SI SI
## 119 2 SI SI
## 120 2 SI SI
## 121 2 SI SI
## 122 2 SI SI
## 125 2 SI SI
## 127 2 SI SI
## 128 2 SI SI
## 129 2 SI SI
## Posee esquema de manejo de plagas y enfermedades
## 3 SI
## 5 SI
## 9 SI
## 11 SI
## 17 SI
## 19 SI
## 20 SI
## 22 SI
## 34 SI
## 35 SI
## 37 SI
## 42 SI
## 63 SI
## 68 SI
## 72 SI
## 80 SI
## 82 SI
## 92 SI
## 101 SI
## 102 SI
## 103 SI
## 107 SI
## 109 SI
## 111 SI
## 113 SI
## 115 SI
## 116 SI
## 117 SI
## 119 SI
## 120 SI
## 121 SI
## 122 SI
## 125 SI
## 127 SI
## 128 SI
## 129 SI
## Ganancias del cultivo de palma Tipo de trabajadores
## 3 2 1
## 5 2 1
## 9 2 1
## 11 2 1
## 17 2 1
## 19 2 1
## 20 2 1
## 22 2 1
## 34 2 1
## 35 2 1
## 37 2 1
## 42 2 1
## 63 2 1
## 68 2 1
## 72 2 1
## 80 2 1
## 82 2 1
## 92 2 1
## 101 2 1
## 102 2 1
## 103 2 1
## 107 2 1
## 109 2 1
## 111 2 1
## 113 2 1
## 115 2 1
## 116 2 1
## 117 2 1
## 119 2 1
## 120 2 1
## 121 2 1
## 122 2 1
## 125 2 1
## 127 2 1
## 128 2 1
## 129 2 1
cat('\n\tCaracteristicas Grupo 1')
##
## Caracteristicas Grupo 1
t(summary(dfsl[[1]]))
##
## Sexo 1:102
## Rango/edad 1: 1
## Sabe leer y escribir 0: 6
## TÃtulo de Educación Superior 0:73
## Nivel de Educación 0: 6
## Categoria de productores - Fedepalma 1:126
## Adopción de tecnologÃa dicotómica 1:111
## Trabaja en actividades diferentes al cultivo de palma NO:94
## Tipo de tenencia 1: 9
## Uso del predio 1:74
## Utiliza sistema de riego NO: 3
## Utiliza fertilizantes NO: 9
## Posee esquema de manejo de plagas y enfermedades NO: 14
## Ganancias del cultivo de palma 1:15
## Tipo de trabajadores 1:84
## grupo Min. :1.000
##
## Sexo 2: 29
## Rango/edad 2:38
## Sabe leer y escribir 1:125
## TÃtulo de Educación Superior 1:21
## Nivel de Educación 1:67
## Categoria de productores - Fedepalma 2: 5
## Adopción de tecnologÃa dicotómica 2: 20
## Trabaja en actividades diferentes al cultivo de palma SI:37
## Tipo de tenencia 2:70
## Uso del predio 2:57
## Utiliza sistema de riego SI:128
## Utiliza fertilizantes SI:122
## Posee esquema de manejo de plagas y enfermedades SI:117
## Ganancias del cultivo de palma 2:68
## Tipo de trabajadores 2:37
## grupo 1st Qu.:1.000
##
## Sexo
## Rango/edad 3:45
## Sabe leer y escribir
## TÃtulo de Educación Superior 2:37
## Nivel de Educación 2:58
## Categoria de productores - Fedepalma
## Adopción de tecnologÃa dicotómica
## Trabaja en actividades diferentes al cultivo de palma
## Tipo de tenencia 3:50
## Uso del predio
## Utiliza sistema de riego
## Utiliza fertilizantes
## Posee esquema de manejo de plagas y enfermedades
## Ganancias del cultivo de palma 3:32
## Tipo de trabajadores 3:10
## grupo Median :2.000
##
## Sexo
## Rango/edad 4:36
## Sabe leer y escribir
## TÃtulo de Educación Superior
## Nivel de Educación
## Categoria de productores - Fedepalma
## Adopción de tecnologÃa dicotómica
## Trabaja en actividades diferentes al cultivo de palma
## Tipo de tenencia 4: 1
## Uso del predio
## Utiliza sistema de riego
## Utiliza fertilizantes
## Posee esquema de manejo de plagas y enfermedades
## Ganancias del cultivo de palma 4:10
## Tipo de trabajadores
## grupo Mean :2.397
##
## Sexo
## Rango/edad 5:11
## Sabe leer y escribir
## TÃtulo de Educación Superior
## Nivel de Educación
## Categoria de productores - Fedepalma
## Adopción de tecnologÃa dicotómica
## Trabaja en actividades diferentes al cultivo de palma
## Tipo de tenencia 5: 1
## Uso del predio
## Utiliza sistema de riego
## Utiliza fertilizantes
## Posee esquema de manejo de plagas y enfermedades
## Ganancias del cultivo de palma 5: 6
## Tipo de trabajadores
## grupo 3rd Qu.:4.000
##
## Sexo
## Rango/edad
## Sabe leer y escribir
## TÃtulo de Educación Superior
## Nivel de Educación
## Categoria de productores - Fedepalma
## Adopción de tecnologÃa dicotómica
## Trabaja en actividades diferentes al cultivo de palma
## Tipo de tenencia
## Uso del predio
## Utiliza sistema de riego
## Utiliza fertilizantes
## Posee esquema de manejo de plagas y enfermedades
## Ganancias del cultivo de palma
## Tipo de trabajadores
## grupo Max. :4.000
cat('\n\tCaracteristicas Grupo 2')
##
## Caracteristicas Grupo 2
t(summary(dfsl[[2]]))
##
## Sexo 1:43 2: 0
## Rango/edad 1: 0 2:10 3:13
## Sabe leer y escribir 0: 2 1:41
## TÃtulo de Educación Superior 0:37 1: 5 2: 1
## Nivel de Educación 0: 2 1:35 2: 6
## Categoria de productores - Fedepalma 1:43 2: 0
## Adopción de tecnologÃa dicotómica 1:41 2: 2
## Trabaja en actividades diferentes al cultivo de palma NO:29 SI:14
## Tipo de tenencia 1: 5 2:34 3: 3
## Uso del predio 1:34 2: 9
## Utiliza sistema de riego NO: 3 SI:40
## Utiliza fertilizantes NO: 6 SI:37
## Posee esquema de manejo de plagas y enfermedades NO:10 SI:33
## Ganancias del cultivo de palma 1: 9 2:15 3:15
## Tipo de trabajadores 1:24 2:12 3: 7
##
## Sexo
## Rango/edad 4:11 5: 9
## Sabe leer y escribir
## TÃtulo de Educación Superior
## Nivel de Educación
## Categoria de productores - Fedepalma
## Adopción de tecnologÃa dicotómica
## Trabaja en actividades diferentes al cultivo de palma
## Tipo de tenencia 4: 1 5: 0
## Uso del predio
## Utiliza sistema de riego
## Utiliza fertilizantes
## Posee esquema de manejo de plagas y enfermedades
## Ganancias del cultivo de palma 4: 3 5: 1
## Tipo de trabajadores
cat('\n\tCaracteristicas Grupo 3')
##
## Caracteristicas Grupo 3
t(summary(dfsl[[3]]))
##
## Sexo 1: 2 2:27
## Rango/edad 1: 0 2:13 3: 7
## Sabe leer y escribir 0: 3 1:26
## TÃtulo de Educación Superior 0:19 1: 9 2: 1
## Nivel de Educación 0: 3 1:16 2:10
## Categoria de productores - Fedepalma 1:29 2: 0
## Adopción de tecnologÃa dicotómica 1:20 2: 9
## Trabaja en actividades diferentes al cultivo de palma NO:21 SI: 8
## Tipo de tenencia 1: 3 2:17 3: 8
## Uso del predio 1:20 2: 9
## Utiliza sistema de riego NO: 0 SI:29
## Utiliza fertilizantes NO: 2 SI:27
## Posee esquema de manejo de plagas y enfermedades NO: 2 SI:27
## Ganancias del cultivo de palma 1: 4 2:10 3:11
## Tipo de trabajadores 1:19 2: 8 3: 2
##
## Sexo
## Rango/edad 4: 9 5: 0
## Sabe leer y escribir
## TÃtulo de Educación Superior
## Nivel de Educación
## Categoria de productores - Fedepalma
## Adopción de tecnologÃa dicotómica
## Trabaja en actividades diferentes al cultivo de palma
## Tipo de tenencia 4: 0 5: 1
## Uso del predio
## Utiliza sistema de riego
## Utiliza fertilizantes
## Posee esquema de manejo de plagas y enfermedades
## Ganancias del cultivo de palma 4: 3 5: 1
## Tipo de trabajadores
cat('\n\tCaracteristicas Grupo 4')
##
## Caracteristicas Grupo 4
t(summary(dfsl[[4]]))
##
## Sexo 1:23 2: 0
## Rango/edad 1:0 2:7 3:8
## Sabe leer y escribir 0: 0 1:23
## TÃtulo de Educación Superior 0: 2 1: 2 2:19
## Nivel de Educación 0: 0 1: 2 2:21
## Categoria de productores - Fedepalma 1:20 2: 3
## Adopción de tecnologÃa dicotómica 1:17 2: 6
## Trabaja en actividades diferentes al cultivo de palma NO: 8 SI:15
## Tipo de tenencia 1: 1 2:19 3: 3
## Uso del predio 1:20 2: 3
## Utiliza sistema de riego NO: 0 SI:23
## Utiliza fertilizantes NO: 1 SI:22
## Posee esquema de manejo de plagas y enfermedades NO: 2 SI:21
## Ganancias del cultivo de palma 1:2 2:7 3:6
## Tipo de trabajadores 1: 5 2:17 3: 1
##
## Sexo
## Rango/edad 4:8 5:0
## Sabe leer y escribir
## TÃtulo de Educación Superior
## Nivel de Educación
## Categoria de productores - Fedepalma
## Adopción de tecnologÃa dicotómica
## Trabaja en actividades diferentes al cultivo de palma
## Tipo de tenencia 4: 0 5: 0
## Uso del predio
## Utiliza sistema de riego
## Utiliza fertilizantes
## Posee esquema de manejo de plagas y enfermedades
## Ganancias del cultivo de palma 4:4 5:4
## Tipo de trabajadores
df_res2 = df_res
cn = c(
'v1' = 'Sex (v1)',
'v2' = 'Range/Age (v2)',
'v3' = 'Can read and write (v3)',
'v4' = 'Higher education degree (v4)',
'v5' = 'Education level (v5)',
'v6' = 'Producer category - Fedepalma (v6)',
'v7' = 'Adoption of Dichotomous Technology (v7)',
'v8' = 'Work in activities other than farming (v8)',
'v9' = 'Type of tenure (v9)',
'v10' = 'Land use (v10)',
'v11' = 'Uses irrigation system (v11)',
'v12' = 'Uses fertilizers (v12)',
'v13' = 'Has pest and disease management scheme (v13)',
'v14' = 'Profits from palm cultivation (v14)',
'v15' = 'Type of workers (v15)'
)
colnames(df_res2) = names(cn)
levels(df_res2$v1) = c('Male','Female')
levels(df_res2$v2) = c('<30','31-50','51-65','66-80','>81')
levels(df_res2$v3) = c('No','Yes')
levels(df_res2$v4) = c('No','Technician or Technologist','Professional')
levels(df_res2$v5) = c('Low','Medium','High')
levels(df_res2$v6) = c('<50 ha','>50 ha')
levels(df_res2$v7) = c('No','Yes')
levels(df_res2$v8) = c('No','Yes')
levels(df_res2$v9) = c('Inheritance succession', 'Owner with title', 'Holder', 'Tenant', 'Co-owner with title')
levels(df_res2$v10) = c('Exc production','Production and residence')
levels(df_res2$v11) = c('No','Yes')
levels(df_res2$v12) = c('No','Yes')
levels(df_res2$v13) = c('No','Yes')
levels(df_res2$v14) = c('0-500000', '500000-1500000', '1500000-3000000', '3000000-6000000', '>6000000')
levels(df_res2$v15) = c('Family, Wage earners', 'Family', 'Wage earners')
colnames(df_res2) = cn
# df_res2 = mutate_all(df_res2, str_to_lower)
head(df_res2)
df_cont = list()
for(i in colnames(df_res2)){
df_cont[[i]] = df_res2 %>%
mutate(v = paste0(i,'_',get(i)),
v2 = paste0(str_extract(i, 'v[0-9]+'),'_',get(i)), # (str_extract(i, '\\((.*)\\)'))
grupo = as.factor(grupo)) %>%
group_by(v, v2, grupo) %>%
summarise(n = n()) %>%
mutate(porc = 100 * n/sum(n))
}
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
## `summarise()` has grouped output by 'v', 'v2'. You can override using the
## `.groups` argument.
df_cont_tot = do.call(rbind, df_cont)
df_cont_tot$ord = nrow(df_cont_tot):1
df_cont_tot
# Los nombres de los grupos son:
# Grupo 1 Población vulnerable
# Grupo 2 Población con liderazgo femenino
# Grupo 3 Población diversificada
# Grupo 4 Población con exclusividad en cultivos de palma de aceite
levels(df_cont_tot$grupo) = c(
'Vulnerable population',
'Population with female leadership',
'Diversified population',
'Population with exclusivity in oil palm crops'
)
graf_porc = df_cont_tot %>%
ggplot()+
aes(porc, fct_reorder(v, ord), fill = grupo)+
geom_bar(stat = 'identity')+
scale_fill_viridis_d()+
labs(x='Percentage',
y='',
fill='Group')+
theme_bw()+
theme(legend.position = 'bottom',
legend.direction = 'vertical',
legend.justification = 'top')
graf_porc
df_cont_ = do.call(rbind, df_cont[1:5])
nom_group = c(
'Vulnerable population',
'Population with\nfemale leadership',
'Diversified population',
'Population with exclusivity\nin oil palm crops'
)
df_radar_tot = df_cont_ %>%
reshape2::dcast('grupo~v2', value.var = 'porc')
df_radar_tot = df_radar_tot[,-1]
df_radar_tot[is.na(df_radar_tot)]=0
df_radar_tot_g = rbind(rep(100, ncol(df_radar_tot)), rep(0,ncol(df_radar_tot)), df_radar_tot)
color = RColorBrewer::brewer.pal(4, 'Set1')
color_alpha = paste0(color, '40')
radarchart(df_radar_tot_g,
cglty = 1, # Grid line type
cglcol = "gray", # Grid line color
pcol = color,
plwd = 2.5, # Width for each line
plty = 1,
pfcol = color_alpha,
vlcex = 1.5)
legend(-1.55, -0.75, horiz = F,
legend = nom_group,
bty = "n", pch = 16,
col = RColorBrewer::brewer.pal(4, 'Set1'),
text.col = "grey25", pt.cex = 3, cex = 1.5)
df_cont_ = do.call(rbind, df_cont[6:10])
df_radar_tot = df_cont_ %>%
reshape2::dcast('grupo~v2', value.var = 'porc')
df_radar_tot = df_radar_tot[,-1]
df_radar_tot[is.na(df_radar_tot)]=0
df_radar_tot_g = rbind(rep(100, ncol(df_radar_tot)), rep(0,ncol(df_radar_tot)), df_radar_tot)
color = RColorBrewer::brewer.pal(4, 'Set1')
color_alpha = paste0(color, '40')
radarchart(df_radar_tot_g,
cglty = 1, # Grid line type
cglcol = "gray", # Grid line color
pcol = color,
plwd = 2.5, # Width for each line
plty = 1,
pfcol = color_alpha,
vlcex = 1.5)
legend(-1.55, -0.75, horiz = F,
legend = nom_group,
bty = "n", pch = 16,
col = RColorBrewer::brewer.pal(4, 'Set1'),
text.col = "grey25", pt.cex = 3, cex = 1.5)
df_cont_ = do.call(rbind, df_cont[11:15])
df_radar_tot = df_cont_ %>%
reshape2::dcast('grupo~v2', value.var = 'porc')
df_radar_tot = df_radar_tot[,-1]
df_radar_tot[is.na(df_radar_tot)]=0
df_radar_tot_g = rbind(rep(100, ncol(df_radar_tot)), rep(0,ncol(df_radar_tot)), df_radar_tot)
color = RColorBrewer::brewer.pal(4, 'Set1')
color_alpha = paste0(color, '40')
radarchart(df_radar_tot_g,
cglty = 1, # Grid line type
cglcol = "gray", # Grid line color
pcol = color,
plwd = 2.5, # Width for each line
plty = 1,
pfcol = color_alpha,
vlcex = 1.5)
legend(-1.55, -0.75, horiz = F,
legend = nom_group,
bty = "n", pch = 16,
col = RColorBrewer::brewer.pal(4, 'Set1'),
text.col = "grey25", pt.cex = 3, cex = 1.5)